Boundary effects on localized structures in spatially extended systems 
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We present a general method of analyzing the influence of finite size and boundary effects on the 
dynamics of localized solutions of non-linear spatially extended systems. The dynamics of localized 
structures in infinite systems involve solvability conditions that require projection onto a Goldstone 
mode. Our method works by extending the solvability conditions to finite sized systems, by incor- 
porating the finite sized modifications of the Goldstone mode and associated nonzero eigenvalue. 
We apply this method to the special case of non-equilibrium domain walls under the influence of 
Dirichlet boundary conditions in a parametrically forced complex Ginzburg Landau equation, where 
we examine exotic nonuniform domain wall motion due to the influence of boundary conditions. 

PACS numbers: 82.40.-e,05.70.Ln 



I. INTRODUCTION 

Models of non-linear spatially extended systems ex- 
hibit a variety of spatial and temporal pattern forming 
phenomena. A subclass of these patterns are spatially lo- 
calized structures [lj that include pulses, solitons, fronts, 
and domain walls. The standard analysis of these local- 
ized structures assumes that, on large length and time 
scales, they can be treated as "coherent objects" p], 
with effective parameters like position, and velocity at- 
tributed to them. A perturbative expansion about this 
isolated coherent object profile is then used to understand 
its response to external forces, interaction with other lo- 
calized structures |2|, |3|, noise, or internal instabilities 
0, . Perturbative calculations reduce the original non- 
linear problem to a series of linear problems that require 
consistency criteria known as solvability conditions for 
their solution. Typically, the solution of a linear equa- 
tion L4> = ipy requires the orthogonality of tp to the zero 
modes ie., (ip,x) = 0, m the null space of the adjoint 
homogeneous problem L'x = 0- 

Often, the symmetries in a particular system are re- 
sponsible for the zero modes of the operators obtained 
after a perturbative expansion. For instance, since a lo- 
calized structure profile and the same profile translated 
infinitesimally are both solutions of the underlying non- 
linear equation, the difference of the two profiles pro- 
vides a zero (neutral or Goldstone) mode. Strictly, the 
zero mode is the derivative of the localized structure pro- 
file, and the underlying symmetry is translation invari- 
ance. Zero modes extracted from symmetry arguments 
may then be employed straightforwardly into solvability 
integrals. 

The argument above, based on translational invari- 
ance, works if the system size is infinite. For a lo- 
calized structure near a system boundary, due to the 
relevant boundary conditions that have to be imposed 
there, the localized structure solution and its infinitcsi- 
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mally translated counterpart are no longer solutions of 
the same equation. Hence, translational invariance is 
broken. Therefore, in this case, one has to contend, not 
only with the incorporation of the boundary data into the 
solvability condition, but also the appropriate treatment 
of broken translation invariance. 

Most treatments of localized structures follow analyti- 
cal techniques that fall in the realm of moving boundary 
approximations 6]. A common feature to these approx- 
imations, for instance, in excitable waves 0, or bistable 
fronts 0, 0, 0| , is the separation of the description of 
the localized structure into an "inner region" and "outer 
region" . The inner region, characterized by short spatial 
scales and fast time scales, captures the internal dynam- 
ics of the localized structure. In contrast, the dynamics 
of the localized structure as a whole is captured by the 
long spatial and time scales comprising the outer prob- 
lem. The solvability integrals in moving boundary type 
approximations occur in the inner problem. Since it is 
the fields in the outer region that mediate the interaction 
with the boundary [TTL fl^ | , the boundary data is not in- 
corporated into solvability conditions arising in the inner 
problem. There are ample situations however, where it 
may not be possible to have separate inner and outer 
regions of a localized structure by manipulating relevant 
system parameters |l4j . In such cases, the boundary data 
must be directly incorporated into the solvability condi- 
tion. 

In this paper, through an appropriately chosen adjoint 
operator L> defined for the semi-infinite system (local- 
ized structure near a boundary), we develop techniques 
that not only include the boundary data into the solv- 
ability condition, but also directly incorporate the effects 
of broken translational invariance into it. We accomplish 
this by extending the definition of the Goldstone mode to 
include the possibility that the corresponding eigenvalue 
be non-zero, with its magnitude dependent on how close 
the localized structure is to the boundary. This leads 
further to a modified solvability criteria. 

As a case study, we develop our techniques in the 
context of reaction-diffusion systems and apply it to 
non-equilibrium domain walls (fronts) found in bistable 
regimes. In bistable reaction-diffusion systems, fronts 
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connecting the two homogeneous steady states can un- 
dergo a bifurcation, called a front bifurcation, where 
a stationary Ising front loses stability to two counter- 
propagating Bloch fronts^). This bifurcation can be 
regarded as an internal instability of the Ising front, the 
localized structure about which a perturbative expansion 
is carried out to obtain the propagating Bloch wall so- 
lution. This bifurcation, also known as the Ising-Bloch 
bifurcation, has been observed in several systems, like 
chemical reactions 0, 0, 0] and also in liquid crystals 

mm. 



In a recent work |12j |. we examined the influence of 
boundaries on Ising-Bloch fronts in a FitzHugh-Nagumo 
(FHN) reaction diffusion model. We were able to de- 
rive order parameter equations (OPE) for front dynam- 
ics, where the fronts were perturbed by the imposition of 
Dirichlet and possibly other boundary conditions at the 
boundaries. This derivation for the two component FHN 
model required restrictive assumptions about the relative 
size of the fronts for the two concentration fields, allow- 
ing for the use of moving boundary approximation like 
singular perturbation methods detailed in 0, [t| . These 
singular perturbation techniques are quite versatile, pre- 
dicting exotic phenomena like front reversal, trapping, 
and oscillation at the boundary. However, as observed 
earlier, we wish to examine the effects of boundary data 
on localized structures, where moving boundary type ap- 
proximations are not applicable, and the explicit incor- 
poration of boundary data in a solvability condition is 
required. 



In the next section we discuss the extension of the 
solvability condition to incorporate boundary data and 
broken translational invariance via the extension of the 
Goldstone mode in a generic system exhibiting a localized 
structure. In Sec. Ill, we describe the modification of the 
slow manifold of a generic Ising-Bloch front due to bound- 
ary effects. In Sec. IV, we apply our method of solv- 
ability condition extension to study the effects of finite 
size and Dirichlet boundary conditions on the dynamics 
of Ising-Bloch fronts in a parametrically forced complex 
Ginzburg Landau equation (CGLE) 0,mi- An impor- 
tant reason behind this choice is its experimental context, 
modeling Ising-Bloch fronts in Liquid crystals subjected 
to rotating magnetic fields m> mi • Liquid crystal sys- 
tems are ideal candidates to study boundary effects, as 
lateral boundary conditions may be imposed in a con- 
trolled manner by appropriate electric fields llSl . An- 
other experimental test bed is presented Ref. | 2f| , in the 
form of coupled non-linear electrical oscillators, where 
the application of boundary conditions requires a minor 
and straightforward variation of the original circuit. In 
Sec. V we discuss in detail the implications of the derived 
order parameter equations for the parametrically forced 
CGLE. In Sec. VI we present our conclusions. 



II. GOLDSTONE MODES AND SOLVABILITY 
CRITERIA 

Consider a general non- linear PDE, 

d t U = CU + N(U), (1) 

where U(x,t) is the solution vector, £ are the linear 
terms, and N(U) are the non-linear terms. Let Uo(x) 
be a stationary localized solution of Eq. QJ, with the 
asymptotic behavior U(x) — > 0;x — > ±oo. In principle, 
Uq(x) also encompasses uniformly translating localized 
structures, which are stationary in a co-moving frame. 

Due to translational invariance in the system, one has 
A(x) = Uo x , the derivative with respect to x of the local- 
ized structure profile, as the zero eigenvalue (neutral or 
Goldstone) mode of the operator £ = £ + N'(Uq). Also, 
it is reasonable to expect that due to translational invari- 
ance £^ has a corresponding zero eigenvector, given by 
the solution of £ jf A' f = 0. A detailed discussion of this 
issue may be found in |13| and the references therein. 

While examining the stability of A = Uq x to perturba- 
tions, which may include a small external perturbation 
p(U,x) added onto Eq. QJ, one obtains, 

[C + N'(U )]6U = f; 
f = d t (SU)-{N"(U )(SU) 2 /2+p(U ,x) 

+ p'(U ,x)5U + p"(U , x)(SU) 2 /2 + 0[(SU) 3 }},(2) 

where SU is the small deviation from the localized struc- 
ture profile. Realizing that the operator £ = £ + N'(Uo) 
has a Goldstone mode, the solvability of Eq. requires, 

CM f )=0. (3) 

The brackets indicate an inner product or the projection 
of the dynamical terms / onto the Goldstone mode (its 
corresponding adjoint) A^ . Equation. © represents the 
generic response of a localized structure to a wide variety 
of perturbations, both internal and external. From an in- 
formal and intuitively appealing point of view, the Gold- 
stone mode with its associated zero eigenvalue is a slow 
(relevant) mode, which coupled with other slow modes 
in the system, should dominate the dynamics. The pro- 
jection in Eq. J2J is a formal prescription to capture this 
slow dynamics. 

Let a localized structure be located near a boundary 
at x = —I, with the origin fixed at the position of the 
localized structure. Although, A* is still a solution of 
£i A^ = in this case, it does not assume the homoge- 
neous boundary value A'(—l) = 0. Consequently, A^ is 
no longer the zero eigenvector of the adjoint homogeneous 
problem in the semi-finite interval [— I, oo]. However, we 
still expect A^ to play a central role in the dynamics 
of the localized structure, all be it in a slightly modi- 
fied form A\ = A^ + 5A\. The subscript I denotes the 
proximity of the localized structure to the boundary, and 
SA] is a proximity dependent correction to A^ . We re- 
quire that in the limit I — ► oo, Aj — > A\ and 8A\ — > 0. 
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This requirement is reasonable on physical grounds. The 
slow dynamics of the localized structure far away from 
the boundary involves as a relevant constituent by 
virtue of it being a slow mode. As the localized structure 
gradually nears the boundary we still expect A\ in its 
modified form Aj, to be the relevant (slow) constituent 
of the dynamics. 

Aj may be determined in two possible ways. Firstly 
we may extract Aj as the solution of 

£^Aj=0, Aj(-l) = 0, A](oo) = 0, (4) 

with the implication that Aj = A^ + SAj is still a zero 
eigenvector in the finite system. Or we may extract A, 
as a solution of 

£U| = A,4, A\(-l)=Q. A»H=0. (5) 

Thus, as the localized structure gradually closes in on a 
boundary the zero eigenvector A'' is modified to A\, and 
the zero eigenvalue gradually migrates away from zero, 
assuming the value A;. Hence, as / — > oo, A; — > 0, and 
4 - A. 

The first scenario is easily discarded using unique- 
ness arguments. If Eq. J3J is obeyed, then 5A\ should 
obey, £HA\ = 0,6Aj(-l) = -A^(-l),SAj(oo) = 0, 
with the unique solution 6 A] = —A^. Therefore, since 
A] = A^ + 5 A], Eq. (@J only has the trivial solution 
A] = (the uniqueness of homogeneous and inhomoge- 
neous problems involving linear differential operators on 
semi-infinite intervals can be proved by a transformation 
that takes the semi-infinite interval into a finite inter- 
val, followed by the utilization of theorems on uniqueness 
available for finite intervals. We provide a proof in Ap- 
pendix A for the CGLE that is studied in detail in later 
sections. Moreover, such a transformation may also be 
applied to operators with an asymptotic structure simi- 
lar to that of the CGLE). This leads us to conclude that 
the modification of A' in a finite system is appropriately 
represented by Eq. (JSJ). 

For arbitrary functions u (not the field U in Eq. Q) 
and i>, and using integration by parts, we have, 

(£u,v) — (u, £^v) + v(b)u x (b) — v{a)u x {a) 

+ v x (a)u(a) - v x (b)u(b), (6) 

where we assume for simplicity that £ is a reaction- 
diffusion type operator comprised of second order differ- 
ential terms only, x = a and x = b are arbitrary bound- 
ary points. If needed, one may evaluate surface terms for 
more general operators using integration by parts. 

For the localized structure a = — I and b = oo. We 
invoke Eq. j£| and substitute v = Aj, u = SUi (the sub- 
script I denotes that 5U is now considered in a finite 
system) in Eq. @, to obtain, 

{£5U h A\) = (f, Aj) = (SUi, XiA\) + Al(-l)6Ut(-l) 
- 4,(00)^(00). (7) 



This is the sought after finite system extension of the 
solvability criteria Eq. (0. Also, as I — > oo, Eq. J2J re- 
duces to (/,At) = 0. Since £ is obtained by linearizing 
about the localized structure Uq(x), SUi(-l) is simply the 
difference U(—l) — Uo(—l), where U(—l) is the Dirich- 
let boundary value imposed on field U, the solution of 
Eq. Q. 

The extension Eq. JJJ), tailored to incorporate non- 
homogeneous Dirichlet boundary conditions on the field 
U, is not unique. For instance, one may consider the ef- 
fects of non-homogeneous Neumann boundary conditions 
on the field U by requiring that Aj obeys 

^4 = ^4, 4(-o = o, 4m = o. (8) 

Here, the derivatives, rather than Aj itself, assume zero 
values at the boundary. Furthermore, an extension Aj for 
a general set of homogeneous boundary conditions, with 
homogeneous Dirichlet and Neumann boundary condi- 
tions as special cases, may also be developed. Next, we 
apply the techniques and criteria developed so far to an- 
alyze non-equilibrium Ising-Bloch fronts, as the fronts 
interact with the system boundary. 

III. BOUNDARY EFFECTS IN A GENERIC 
ISING-BLOCH SYSTEM 

Ising-Bloch fronts provide an interesting arena to apply 
the methods developed in the last section. Along with 
the usual Goldstone mode associated with translational 
invariance, the slow manifold for Ising-Bloch fronts also 
includes a spatially localized slow mode responsible for 
the Ising-Bloch bifurcation ^.pll 12^ . Chirality preserv- 
ing stationary Ising fronts |14j . bifurcate into a pair of 
chirality broken, counter-propagating Bloch fronts. The 
slow manifold for Ising-Bloch fronts comprised of the 
Goldstone and chirality breaking modes, manifests itself 
in the form of order parameter equations (OPE) 0,0, 
for the order parameters, front velocity and front posi- 
tion. The front velocity is a measure of the effects of 
the chirality breaking mode. The Goldstone mode cap- 
tures front translations by infinitesimal changes in the 
front position, the other order parameter. We seek the 
coupling between these order parameters induced by the 
boundary data and broken translational invariance. 

A generic Ising front denoted by Uq(x), gives the Gold- 
stone mode Uq x . Close to the Ising-Bloch bifurcation 
threshold, propagating Bloch wall solutions are regarded 
as perturbations of the stationary Ising wall solution [J^ . 
The front velocity c controls the strength of these per- 
turbations. Therefore, expanding the deviation SU in 
powers of c, we have, 

U b = U + 5U 

= U +cSU 1 + c 2 5U 2 + c 3 6U 3 + .., (9) 

with the perturbed Bloch wall solution Ub- 
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For convenience we transform into a frame of reference 
moving along with the Bloch wall. This transformation 
amounts to d t {SU) — > d t (SU) — c(Uq x + SU X ). Invok- 
ing Eq. @ and substituting into it the expansion of SU, 
while at the same time disregarding the influence of any 
external perturbation p(U, x), we obtain, 

£[c5Ux + c 2 SU 2 + c 3 SU 3 ] = d t {c5U{) 
- c[U 0x + cSU lx + c 2 8U 2x ] - c 2 N 2 - c 3 N 3 + ■ ■ (10) 

N 2 and A3 represent the coefficients of second order and 
third order velocity terms respectively. 

Equating terms which are first order in velocity c in 
Eq. (jTU|) . we obtain, 

£SUi + U ax = 0. (11) 

This means that £ has a double zero eigenvalue at 
the Ising-Bloch bifurcation threshold 0,112 ■ Therefore, 
along with the zero Goldstone mode, we have another 
eigenvalue that passes through zero at the bifurcation. 
The Goldstone mode Uo x and the generalized eigenvec- 
tor 5Ui obtained from Eq. 1)11(1 . span the slow manifold. 
The chirality breaking mode is then constructed as a lin- 
ear combination of these two modes [5j. 

Employing the projection criteria Eq. © for an Ising- 
Bloch front close to the bifurcation threshold, ie., the 
solvability of Eq. (|10|) , results in, 

(SU u A^d t c = c(U 0x> A^ + c 2 (SU lx + N 2 ,A^) 

+ c 3 {5U 2x +N 3l A^) + --- (12) 

This is the generic OPE for the velocity of Ising-Bloch 
fronts close to the bifurcation threshold. The particular 
form of the inner products in Eq. (|12|) is system specific. 
If one assumes further symmetries in the system, for ex- 
ample U — ► — U, inner products that are coefficients of 
even powers of the velocity in Eq. I|12l) vanish, resulting 
in the normal form of a pitchfork bifurcation. The inner 
product (Uo x ,A' i ) in Eq. i|12|) controls the distance from 
the Ising-Bloch bifurcation threshold, where for consis- 
tency (Ising-Bloch bifurcation is a pitchfork) it is further 
required that (Uq Xi A^) ~ c 2 , dtc ~ c 3 UH^]. Hence, all 
the terms in Eq. (|12|) are of size c 3 . 

We invoke the extended solvability criteria Eq. (J7J to 
evaluate the effects of boundary data on the dynamics of 
Ising-Bloch fronts. For generic Ising-Bloch fronts inter- 
acting with boundaries where Dirichlet data is present, 
the extended solvability criteria assumes the form, 

{5U lu A\)d t c = c(U 0x ,A])+c 2 (6U llx +N 2 ,Al) 

+ c 3 {5U 2lx +N 3 ,A\) 

+ Xi(cSU u + c 2 6U 2t + c 3 6U 3l + ■ ■ • , A\) 

+ A\ x {-l)6Ui{-l) - 4 x (oo)^,(oo).(13) 

In contrast to earlier works 0,0, E2 focused on the effects 
of external perturbations, p(U,x), on the slow manifold, 
the constituent modes of the slow manifold require appro- 
priate modifications in order to capture the effects arising 



due to confinement by boundaries. While, the modifica- 
tion of the adjoint Goldstone mode A^ to A\ is generic to 
any confined localized structure, or alternatively, a local- 
ized structure being considered in the vicinity of system 
boundaries, the modification of the generalized eigenvec- 
tor SUi to 5Uu is a unique characteristic of Ising-Bloch 
fronts. 

Simplifications to the slow manifold Eq. (|13[) arc made 
by the following observations. Consider the term, /o = 
Xi(cSUii+c 2 SU 2 i+c 3 SU 3 i + - ■ ■ , A}), on the right hand side 

of Eq. (T3JI. The inner product /1 = Xi(c5Uu, A\) has the 
largest contribution since it involves the first power of the 
velocity c. Now, as mentioned before, all terms should 
be of size c 3 , a requirement imposed for the Ising-Bloch 
bifurcation to be a pitchfork. Therefore, fi ~ A;c ~ c 3 , 
implying A; ~ c 2 . Moreover, the size of A; is controlled 
by the distance of the Bloch fronts from the boundary. 
If the front is far away from the boundary, that is, if 
A; ~ 0(c 3 ), then /1 ~ C(c 4 ), and its contribution to 
Eq. (|13fl can be neglected. As the front moves towards 
the boundary, so that A; ~ c 2 , then f\ ~ c 3 contributes 
to Eq. (|13l) . and the ensuing front dynamics. If the front 
gets too close to the boundary, ie., A; ~ c, then fx ~ c 2 , 
and the scaling requiring that all the terms be of size c 3 
breaks down. In other words, if A; ~ c, the effects of the 
boundary are too strong for them to be accurately con- 
sidered as small perturbations on the dynamics of Ising- 
Bloch fronts. Consequently, the size of A; serves as a 
measure of the strength of the boundary perturbation. 
In light of the present discussion, Eq. (|13|) simplifies to 

{5U xh A\)d t c = c(U 0x ,Al) + c 2 (SU Ux + N 2 ,Al) 

+ c 3 (8U 2lx +N 3 ,A\) 

+ Xi(cSU u ,Aj) 

+ Al(-l)SUi(-l). (14) 

The surface terms at infinity contribute zero, since by 
construction A', (00) = 0. 

IV. BOUNDARY EFFECTS IN THE 
PARAMETRIC ALLY FORCED CGLE 

The CGLE reads, 

d T F = (7 + iv)F - \F\ 2 F + fiF* + d 2 x F + a. (15) 

Equation. i|15|) and its generalizations US E-l have been 
thoroughly analyzed in the context of the Ising-Bloch bi- 
furcation. The field F may be regarded as the amplitude 
of diffusively coupled auto oscillators that oscillate above 
the Hopf bifurcation threshold determined by the param- 
eter 7. fi represents the strength of parametric forcing at 
twice the natural frequency, and v is the detuning. The 
parameter a, which models forcing at the natural fre- 
quency of the system, breaks the [F — > —F) symmetry. 
As a result, the pitchfork normal form of the Ising-Bloch 
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bifurcation for a — unfolds into a saddle node for a 
non-zero a. 

We briefly recount the results of p| concerning the dy- 
namics of Ising-Bloch fronts in the parametrically forced 
CGLE valid for an infinite system. This lays down the 
framework for the subsequent consideration of finite sys- 
tem sizes and boundary effects. 

For a — and in the bistable regime deter- 
mined by the constraints, \v\ < fx, 7 > —\J fx 2 — v 2 , 
Eq. (|15|l possesses a stationary Ising wall solution Fj = 
Y / Ktanh(-\/K/2X)e 1 ^. Here k = 7 + \/ [i 2 — v 2 and is 
obtained by solving sin(2(/>) = v j [i. Bloch wall solutions 
of Eq. I|15|l are then obtained as a perturbation to the 
Ising wall, 



Fb(x, t) = y / K[tanh(x) + u(x, t) + iw(x, t)] 



where the space-time scaling t = kt/2, x 
introduced by the authors, resulting in, 



Substituting these vectors into Equation. (|12|l gives 




£ 



D 1 -Av/k 
L» 2 -3 + 47/k 

£>! = <9 2 + 2-6tanh 2 (.T), 
£ , 2 = ^ + l-2tanh 2 (a;), 





" 3u 2 + w 2 




u 3 + uw 2 


2 tanh(cc) 




-2 






2uw 




w 3 + wu 2 



For clarity and continuation of the conventions used in 
the previous sections, we stress the following points. 
Firstly, we recognize that SU = {u,w} T . Secondly, SU 
obeys 



d t SU = £SU + N, 



(17) 



which when compared with Eq. leads to the realiza- 
tion that N = N"(U a )(SU) 2 /2 + 0[{SU) 3 . Thirdly, £ 
is obtained by linearizing about the solution Uq(x). In 
the present case the stationary solution is the Ising wall 
Fi(x) = y/n tanh(x)e^, and Uq(x) = tanh(x), where the 
constant factor \fke % ^ should be dropped if the pertur- 
bation SU = {u,w} T is defined through Eq. (|Td|) . 

For the specific case of the parametrically forced 
CGLE, one has @, 



SUi 



and 



^/u(x)- Ji 2 (a:) 
l^sech(x) 





sech 2 (x) 












(18) 



i^^sech 2 (a;) 
sech(x) 



27(/x c - fj,)fi c 
d * c = ^2 c 
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+ 0.36 c 3 . (19) 



Eq. ((TU|) possesses three stationary states, two counter- 
propagating Bloch walls and a stationary Ising wall. 
These steady states exchange stability via the Ising- 
Bloch bifurcation at the critical bifurcation parameter 
3/i c = \/§v 2 + 7 2 . The components of the vectors 
5U = cSUi +c 2 5U2 + ■■, Uq and A' , in an infinite system, 
exponentially decay to zero as one moves away from the 
front both to the left and to the right. This signifies that 
Ising and Bloch walls are localized structures that are 
not influenced by boundary conditions imposed on either 
boundary sufficiently far away. Furthermore, no explicit 
dependence on x in Eq. (|19f) indicates translational in- 
variance, a residue of infinite system size. 

We now calculate A\ and the associated value of A; . A\ 
satisfies the boundary conditions A\{— I) = 0, Aj(oo) = 
(homogeneous problem), since we wish to examine the 
influence of Dirichlct boundary conditions on U (non- 
homogeneous problem). Close to the bifurcation thresh- 
old determined by the magnitude of /i c — /i, the operator 
£ jf has the form 



Di 

-I//7 D 2 
£\ + (fi-n c )£l 



27^ c (/x - n c ) 
4 7 2 





»h - 1 



(20) 



The operator £\ is a perturbative correction to the op- 
erator £\, since fi — [i c ~ c 2 . Hence, we first examine £\ 
the dominant term in £'. 

The operators D\ and D2 populate the diagonals of £[, 
and possess zero eigenvectors given by Z\ — sech 2 (x) and 
Z2 = sech(a;) respectively, in an infinite system. These 
eigenvectors satisfy the constraint of being zero at posi- 
tive and negative infinity. Imagine a traveling Bloch front 
sufficiently distant from the boundary, where Dirichlct 
boundary conditions are imposed. The front does not 
sense the boundary and the condition D\Z\ = D2Z2 = 
holds. This is because the solutions Z\ and Z2 exponen- 
tially approach zero on either side of the front. As the 
front closes in on the boundary, such that it is barely 
able to sense it (Z\ and Z2 have small finite values at 
the boundary), the eigenvectors Z\ and Z2 are modified 
to Z\i and Z21 by constraining them to have zero values 
at the boundary. Meanwhile, in a semi-infinite or finite 
domain, the only solutions to D\Zu — D2Z21 = which 
have a zero value at both boundaries are the trivial so- 
lutions Zu = Z21 = (uniqueness arguments). Hence, 
requiring that the solutions Zu (Z21 ) are only slight mod- 
ifications of Zi(Z2) and are not trivial zero solutions de- 
mands that these solutions obey D\Zu = XuZu and 

D2Z2I = A2/^2/- 

Figure. ^a) shows the plot of Z\ in grey, where the left 
boundary is at a finite distance / from the peak. Z\ has a 
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finite nonzero value at the boundary. We require that the 
modified eigenvector Zu have a zero value at the bound- 
ary and not be all that different from Z\ elsewhere. We 
make the ansatz that this can be accomplished by sub- 
tracting from Z\ its image to the left of the boundary. 
Therefore, we have, Z\i — sech 2 (x) — sech 2 (x + 21). Fig- 
ure, mb) shows a good agreement between our guess and 
the actual numerically evaluated Zu. This is so because 
in the asymptotic limit exp2x >> 1, D\ — <9 2 — 4, and 
the image is approximately a zero eigenvector of this op- 
erator in the same limit. 

Introducing images into a semi-infinite problem is by 
no means a coincidence. Images are a common occur- 
rence whenever boundary data is involved. For the ex- 
tension A] (correspondingly Zu and Z%) to assume a 
zero value at the boundary, the introduction of the im- 
age becomes a natural necessity. Furthermore, we wish 
to stress that the concept of images is quite general in 
its utility. Extensions of Goldstone modes can be readily 
obtained for other systems, with linear operators having 
similar properties of exponential decay asymptotics. 

An upper bound, A{,, on the eigenvalue Xu, is easily 
obtained by a variational principle, given by, 

| An | < |Al,| = {Z l uD l Z ll )/{Z lu Z ll ). (21) 

A more refined variational guess of Zu may be made 
by introducing an extra parameter a±. Consequently, 
we have Zu = exp (aia:)[sech 2 (a;) — sech 2 (x -I- 2/)]. Ma- 
nipulation of this parameter provides a better guess of 
the change in shape of the peak in the actual modi- 
fied eigenvector Zu- Figure. Ufa) compares the numeri- 
cal and variationally calculated eigenvalues as a function 
of the distance / of the front from the boundary. The 
dashed curve represents the numerically calculated eigen- 
values of D\. The thin curve depicts the variationally 
calculated eigenvalues with Zu = sech 2 (a;) — sech 2 (a; + 
21). The squares signify a better variational calcula- 
tion of the eigenvalues using Zu = cxp (aia;)[sech 2 (a;) — 
sech 2 (x -I- 2/)]. An improved guess of Z 2 i, and eigen- 
value X21 for the operator D2, similarly involves taking 
Z12 = exp (a2cc)[sech(a;) — sech(a; + 21)]. Depicted in 
Fig. OJb) are the eigenvalues X21 , numerically calculated 
(dashed curve), variationally calculated with respective 
guesses Z21 = sech(x) — sech(x + 21) (thin line), and 
Z21 = exp (a,2x)[sech(x) — sech(x + 21)] (squares). 

The numerical calculation of the eigenvalues Xu and 
A2Z involved using a standard QR algorithm on the ma- 
trix obtained by a finite difference approximation to the 
operators D\ and Z?2- The grid spacing was adjusted 
until we obtained convergence. The eigenvectors were 
calculated using inverse iterations, with the number of 
iterations optimized for convergence. 

The first row in the matrix representation of the ad- 
joint operator Eq. I|20|) consists only of the operator D\. 
Therefore, since &A\ = XiA}, we immediately obtain 
A; = Ai;. We recall that in the limit of infinite front 
distance from the boundary I — ► 00, we have A\ — > A^ . 




FIG. 1: (a) Shows the plot of Z\. The peak is at a dis- 
tance of I = 2 from the boundary, (b) The squares rep- 
resent the numerically obtained Zu. The analytical guess 
Zu = sech 2 (:r) — sech 2 (a; + 21) is the solid line. 



Combining this asymptotic limit constraint with the re- 
quirement that the sought after eigenvector has zero val- 
ues at both boundaries, we obtain, 



(22) 

A more rigorous derivation involving a step by step con- 
sideration of the operators L[ and L 2 in a perturbative 
scheme also yields Eq. (|22f) . 

We now focus on incorporating the effects of the Dirich- 
let boundary values Xb and Yb, the values of the real 
and imaginary components of the field F in Eq. Q15p. 
into the dynamics of fronts close to the boundary. Bloch 
walls are perturbed Ising walls, with the perturbation 
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FIG. 2: (a) Comparison of variational and numerical calcu- 
lations of Aii (b) Similar comparison of A2; calculated using 
numerical and variational techniques. 



SUi. The boundary value of this perturbation SUi(-l) 
is obtained by fixing F(—l) — X b + iY b and subtracting 
from it the value that the Ising wall assumes Fj(—l) = 
^Ktanh(-Z)e^. Recalling Eq. JTHJ), and SU = {u,w} T , 
we obtain, 



(X b cos(<j)) + Y b am(<f>))/y/K+tanh(!) 
(y 6 cos(0) - X b sm(cf))) / y/U 



(23) 



OPE 



To extract a reduced description of the influence of 
Dirichlet boundary conditions on the motion of Ising- 
Bloch fronts, we invoke Eq. Q14|). and substitute into it 



the explicit forms of A\ and A; derived in the previous 
section. Consider the term fx = Xi(c8Uu, Aj) on the 
right hand side (RHS) of Eq. ([Q. For the CGLE, as seen 
in Eq. (1221) , the first component of A\ , denoted by, A' n , is 

smaller by a factor of c 2 than the second component A l2 . 
This is so because [i c — fi ~ c 2 . Hence, while evaluating 
fi , we need only consider the inner product of the second 
component of the generalized eigenvector, 6Uu, denoted 
by SU112, and A\ 2 . The generalized eigenvector 5U\ is 
known Eq. (|19|l . and its finite system modification SUu 
needs to be evaluated (only the second component SU112) 
to evaluate the inner product in /j . 

To evaluate 5U\i2 we recall that Z 2 = sech(a;), with 
D2Z2 = 0. The second component of 8U1, is given by 
SU12 — [87/97rz/]sech(a;). Hence, D2SU12 = 0. In a con- 
fined system with the left boundary at x = —I, Z2 is 
modified to Z21 — sech(a;) — sech(x + 21), requiring that 
the homogeneous boundary condition, Z%i{—1) = 0, holds 
good. In the confined system 5Ui 2 is modified to 6Uu2- 
However, to obtain SUm, the requirement that it obeys 
the inhomogeneous boundary condition cSUu2{—l) = 
5Ut 2 (—l), since SUi = c5Ui + 0(c 2 ), needs to be imposed. 
Therefore we construct SUu2(x) — cSUi 2 — f3sech(x + 21), 
followed by imposing the inhomogeneous boundary con- 
dition c6Uu2(—l) = 5Ui2{—l), to evaluate (3. After doing 
so, we have, 



C87 

cSUu2 = n Z21 



SUi2(-l) 
sech(Z) 



sech(x + 2Z). (24) 



We, finally have the ingredients to calculate all the 
inner products in Eq. (|14p. The bulk of the bound- 
ary influence, we contend, is captured by the interplay 
of the terms, c(Uq x , Aj), \i(c5Uu, A\), and the surface 
term A\ x {-l)8Ui{-l) in Eq. JSJ. Therefore, although, 
strictly speaking, the inner products containing higher 
order terms c 2 {SU Ux + N 2 ,Aj), and c 3 {SU 2 i x +N 3 ,Aj), 
in Eq. (|14() . should be evaluated in the finite domain 
[— I, 00], we approximate them by taking the inner prod- 
uct in the infinite interval [—00,00]. 

Performing all the inner products in Eq. I|14|l and re- 
arranging the terms, we obtain 



d t c = 



27(fi c - fi)fi c 



4j 2 



c + A(C 



" 87 " 








9nv 





9ttv 



I67 
81(/x 



tanh(l)secb.(l)6Ui 2 (-l) 



4 7 2 



tanh(Z)*f7u(-0 



I67 



2fcosech(2Z). 



(25) 



In deriving Eq. (|25l) we have used Z\i = sech 2 (x) — 
sech 2 (x + 21) and Z21 — scch(a;) — scch(a; + 21), where 
A ; = X u is given by Eq. and p = 0.36 Eq. ifT^I . 

Equation. (|25|l along with dtl — —c represents the cou- 
pling of the two degrees of freedom, front velocity c and 
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position I, by the influence of Dirichlet boundary condi- 
tions imposed at the boundary. As required, in the limit 
of infinite front distance from the boundary Eq. I|25|) re- 
duces to Eq. 

We now examine the consequences of the coupling of 
the front velocity and position close to the boundary. 
Firstly, we report the findings of our numerical simula- 
tions of Eq. (|15J) . which is a system with infinite degrees 
of freedom. Secondly, we corroborate these findings by 
solving the reduced, two degree of freedom OPE we have 
derived. 

We performed numerical simulations of Eq. (|15|l , where 
Bloch fronts were created at infinity (far from the bound- 
aries) and launched towards a boundary. The velocity of 
these Bloch fronts was chosen to be one of the steady 
states of Eq. (|19l) resulting in uniform front translation 
with this velocity until the fronts closed in on the bound- 
ary. Near the boundary, contingent upon the Dirich- 
let boundary value imposed, the incoming Bloch fronts 
were either trapped or bounced back. Bloch fronts that 
bounce evolve into the counter-propagating Bloch front 
near the boundary and move away. Trapped Bloch fronts, 
as opposed to bouncing Bloch fronts, evolve into non- 
trivial steady state solutions (See Ref.jl^l) of the CGLE 
Eq. dSJ. 

We summarize our numerical observations of Bloch 
front behavior as a function of the boundary conditions 
Xf, and Yb in Figure. 0. This phase diagram in the plane 
of boundary values reveals a curve separating regions of 
bouncing and trapped fronts represented by diamonds. 
We compare these results with the transition curve pre- 
dicted by the reduced model Eq. (|25|) . plotted as the 
dashed curve in Figure. |3| The plots show a good agree- 
ment (within 0.5%) between the two transition curves. 
This is a striking result considering the fact that in calcu- 
lating A\ and A; we have employed approximate vectors 
Z u and Z 2 i- 

Bouncing fronts gradually slow down as they near the 
boundary, attain zero velocity at a certain critical dis- 
tance from it, and finally move away as the sign of the 
velocity flips. As we change the boundary values and 
get closer to the transition curve, bouncing fronts attain 
zero velocity at a much smaller critical distance from the 
boundary, until eventually right at the transition curve 
they reach the point of closest approach to the bound- 
ary. As we cross the transition curve and move into the 
trapping region, approaching fronts no longer attain zero 
velocity close to the boundary, their velocity never flips 
sign, and hence they never bounce. The distance from 
the boundary of the point of closest approach depends on 
where exactly on the phase diagram the transition curve 
is crossed. 

The agreement between the transition curves obtained 
from the full model Eq. ljTo|) and the reduced model 
Eq. 125|) is better when the point of closest approach is 
further away from the boundary. This is because, as de- 
tailed earlier, the vectors Zu and Z21 are better approx- 
imations to the actual solutions of D\Zu = \\Z\\ and 




-0.6 -0.55 -0.5 -0.45 -0.4 -0.35 -0.3 

yb 

FIG. 3: The transition curve for the full model Eq. 1)15^ plot- 
ted using squares, the same curve obtained from the reduced 
OPE Eq. 1251 . plotted as a dashed line. Here, v — 0.3, 
7 = 1.0, /i = 0.448. 



D2Z21 = \2Z2u further away from the boundary. Con- 
sequently, a better guess of these vectors, valid close to 
the boundary, should improve the agreement between the 
transition curves, even if, the point of closest approach is 
closer to the boundary. However, the approximate vec- 
tors we use are sufficient for the purpose of establishing 
the usefulness of our general method that accounts for 
the broken translational invariance in a spatially finite 
system through the extension of solvability conditions. 
Our method incorporates into it the eigenvalue A;, the 
most direct measure of broken translational invariance, 
which can be obtained accurately via a variational prin- 
ciple using relatively crude guesses for the eigenvectors. 

We now, by examining Eq. (|25[l in more detail, ex- 
tract the mechanism behind the transition from bounc- 
ing to trapped fronts as Dirichlet boundary conditions 
are changed. Figure, ^a) shows the nullclines, invariant 
manifold, and trajectories of Eq. (|25|l inside the bounc- 
ing region of the phase diagram. A saddle, present at the 
point of intersection of the nullclines, controls the flows in 
this bouncing regime. Far away from the boundary, sit- 
uated at x = in the plot, the nullclines are three paral- 
lel straight lines that represent two counter-propagating 
Bloch wall steady state solutions, and a stationary Ising 
wall solution of Eq. i|19fl . The bouncing involves the 
Bloch front initially flowing towards the saddle. There- 
upon, influenced by the unstable manifold, the front flows 
away. 

Figure. 01b) still depicts flows inside the bouncing re- 
gion, but much closer to the transition curve. In this 
regime bouncing and trapped fronts can coexist. The in- 
variant manifolds demarcate two basins, one of attraction 
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FIG. 4: (a) The plot deep inside the bouncing region, the 
nullclines are thin black curves, the thick curves correspond 
to the trajectories in the phase plane, and the invariant man- 
ifolds are plotted as dashed lines. Here, v = 0.3, 7 = 1.0, 
H = 0.448, X b = -1.116, and Y b = -0.4262. (b) Plot still in 
the bouncing region, but close to the transition curve. The 
same plotting scheme and parameters used, with boundary 
values X b = -1.112, Y b = -0.4262. 



FIG. 5: (a) Plot in the trapping region close to the transition 
curve. The same plotting scheme and parameters used, with 
boundary values X b = -1.11, Y b = -0.4262. (b) The plot 
deep inside the trapping region, the nullclines are thin black 
curves, the trajectory is the thick curve. Here, v = 0.3, 7 = 
1.0, fi = 0.448, X b = -1.09, and Y b = -0.4262. 



towards the boundary, and the other of repulsion away 
from it. Inside the repulsion basin all incoming Bloch 
fronts bounce with the same mechanism as in Fig. QIa) . 
All the flows in the attraction basin are directed towards 
the system boundary, with no possibility of a bounce. 
Figure, ^b) shows both bouncing and trapped Bloch 
front trajectories in their respective basins. We reported 
on the the coexistence region in our numerical study of 
Eq. (|15fl in Ref . jl^] • Here, we have provided an analytical 
explanation of this phenomena. 

The flows in the trapping region close to the transi- 
tion curve are shown in Figure. EJa) ■ Trapped Bloch 
fronts, created at infinity and on the upper branch of the 



nullcline (corresponding to one of the steady states of 
Eq. (|19|l % ). lie inside the basin of attraction towards the 
boundary. Consequently, the transition from bouncing to 
trapped fronts is marked by the initial front velocity and 
position moving from the basin of repulsion (Fig. 01b)) 
to the basin of attraction (Figure. EJa)) as the bound- 
ary values are varied. Deep inside the trapping region 
the saddle no longer exists, and we have a sink instead 
(Fig-GJb)). All incoming Bloch front trajectories end up 
at this sink. 

Summarizing, the nonuniform motion of Bloch fronts 
close to the boundary is governed by the fixed point of 
Eq. (|25[) . giving rise to bouncing, trapping, and coexis- 
tence of the two. Well inside the bouncing region this 
fixed point is a saddle. Deep into the trapping region the 
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fixed point changes into a sink. 

VI. CONCLUSION 

We have developed a general method of analyzing the 
influence of broken translational invariance due to finite 
size and boundary effects on the dynamics of localized so- 
lutions of generic non-linear spatially extended systems. 
We apply our method to the special case of a bistable 
reaction-diffusion system, where the localized solutions 
are fronts Eq. (|25[) . The implementation of this method 
involves the extension of the infinite system size limit 
solvability conditions, used to extract a reduced descrip- 
tion of the infinite dimensional system, into solvability 
conditions that account for finite system size and bound- 
ary effects. The extended solvability criteria works by 
naturally incorporating into it the concept of images. As 
a result, the method affords a direct grasp of the broken 
translational invariance in a confined system through the 
calculation of relevant eigenvalues. 

In the special case of Dirichlet boundary conditions im- 
posed on the CGLE, we were able to provide mechanisms 
for Bloch front trapping, bouncing and coexistence of the 
two at the boundary. This nonuniform front motion is a 
result of the coupling of the two degrees of freedom, front 
velocity and position, by the influence of boundary con- 
ditions. We have explicitly derived this coupling by using 
our method of solvability condition extension. The role 
of other types of boundary conditions, either Neumann or 
mixed can be explored in a similar fashion by construct- 
ing a suitable extension of the modified Goldstone mode. 
For example, exploring Neumann boundary conditions 
requires the extension to always have zero derivatives at 
the boundary. This can be accomplished in the CGLE 
or other systems by adding, rather than subtracting, the 
image. 

Finally, we comment on the generality of solvability 
condition extension. In any system, whenever it is possi- 
ble to derive reduced dynamical equations through pro- 



jections on the Goldstone mode, our method can be ap- 
plied to obtain the finite size and boundary effects in 
terms of the modifications of these reduced dynamical 
equations. 
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APPENDIX A 

For the CGLE, consider the operator D\ Eq. (|17fl in a 
semi-infinite interval [— Z,oo]. Using the transformation 
t = 1 — e~^ l+x \ the problem 

D X Y = [<92 + 2-6tanh 2 (.T)]y = 0, 
Y(-l) = ; y(oo) = 0, (Al) 

is transformed to 

_ ft 2 + 6tanh^ + ln(l-t)) 

(1-*) (1-*) 
F(0) = ; Y(l) = 0. (A2) 

Equation. I|A2(1 has a regular singular point at t = 1, 
and thus has a unique solution. Similar considerations 
apply to the operator Z?2- Therefore, homogeneous or 
inhomogeneous problems involving the operator £, which 
is comprised of the operators Dy, and -D2, should have 
unique solutions in a semi-infinite domain. For operators 
that possess exponential decay asymptotics (true for a 
wide variety of models of physically occurring localized 
structures), a transformation of the type used here, can 
always be found in order to prove the uniqueness. 
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